TUIMP-TH-96/77 



Yukawa Corrections to Top Quark Production at 
the LHC in Two-Higgs-Doublet Models 



Hong-Yi Zhou a ' b , Chong-Sheng Li ac and Yu-Ping Kuang a ' b 

CCAST (World Laboratory), P. 0. Box 8730, Beijing 100080, P.R. China, 
b Institute of Modern Physics and Department of Physics, 
Tsinghua University, Beijing 100084, P.R. China * 
c Physics Department, Peking University, Beijing 100871, P.R. China* 
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I. Introduction 



Recently, the existance of the top quark is confirmed by the CDF and DO Col- 
laborations with the mass and production cross-section mt = 176±8(stat)±10(syst) 
GeV a = 6.8t||p&, and m t = 199±H(stat) ± 22(syst) GeV a = 6.4 ± 2.2pb, 
respectively Jl| . This is a remarkable success of the Standard Model (SM) since 
this measured mass is close to the central value predicted by the best fit of the SM 
to the LEP precision electroweak data. There are a number of new and interesting 
issues related to the top quark, such as the precision measurement of the mass, 
width and Yukawa couplings through its direct production and subsequent decay at 
both hadron and e + e~ colliders. At the future multi-TeV proton colliders such as 
the CERN Large Hadron Collider (LHC), ti production will be enormously larger 
than the Tevatron rates, and the accuracy with which the top mass can be measured 
in proton colliders is around 2 GeV||. This makes the determination of production 
cross section at LHC more accurately than at Fermilab Tevatron. Processes related 
to the top quark (the heaviest fermion so far found) may be sensitive to possible 
new physics beyond the SM, especially the non-standard Higgs sector, due to the 
large Yukawa coupling. 

The simplest extension of the SM Higgs sector is the two-Higgs-doublet models 
(2HDM) which can give an enhancement of the Yukawa couplings and may result 
in a considerable correction to the top quark production cross section. Yukawa 
correction to the top quark production in general 2HDM and in the minimal su- 
persymmetric model(MSSM) at Fermilab Tevatron has been carried out in Ref. ||, 
where large corrections are found. Supersymmetric electroweak correction []5| and 
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QCD correction^ to top quark production at the Tevatron are also found to be 
large. At LHC the main production mechanism of top quark is the gluon-gluon 
fusion process gg — > ti. The full electroweak corrections to this process as well as 
qq — * ti in the Standard Model have been calculated in Ref. jQ. In this paper 
we investigate the Yukawa correction to the top quark production by the process 
gg — > ti in a general 2HDM and in the minimal supersymmetric model at LHC. Our 
calculation of the corrections to the S-matrix elements will be given in Sec. II. In Sec 
III, we present our numerical results and discussions of the cross-sections in the SM 
, 2HDM and MSSM . 

II. Calculation of the S-matrix Elements 

The tree-level Feynman diagrams and the relevant Yukawa corrections to gg — > 
ti are shown in Fig.l (u-channel diagrams of Fig. lb and Fig. Id-Fig. lh are not ex- 
plicitly shown) in which the dashed lines in Fig.lc and Fig.le-Fig.lh represent H°, 
h, A, G°, H + and G + (G° and G + are, respectively, Goldstone bosons of Z and 
W + ), while those in Fig. Id represent only H° and h. The relavant Feynman rules 
can be found in Ref.flTj. 

At tree level, the S-matrix element is composed of three different production 
channels(s-,t-, u-channel) as follows: 



M s , = -zg 2 s zf abc T^u(p 2 M Pl )/s (1) 

M* = -^(T 6 T%e(p 4 )^(p 3 )^(p2)7M(^-^4 + m07^(pi)/(t-^) (2) 

M« = M*(p 3 ^P4, T a ^T b , t-«), (3) 
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where T M is given in the Appendix. Instead of calculating the square of the am- 
plitudes explicitly, we calculate the amplitudes numerically by using the method of 
Ref. || . This method greatly simplifies our calculations. 

The 0{am 2 /m 2 v ) Yukawa corrections to gg — > ti are shown in Fig.lc-Fig.lh, 
where Fig. Id alone is gauge invariant and the sum of Fig.lc,Fig.le-Fig.lh is gauge 
invariant. In our calculation, we use dimensional regularization to regulate the 
ultaviolet divergences and adopt the on-mass-shell renormalization scheme. 

We only give the explicit results of the s- and t-channel contributions to the 
total Yukawa correction. The u-channel results can be obtained by the following 
substitutions: 

Ps <-> Pa, T a «-> T fe , i u. (4) 
Fig.lc leads to the s-channel vertex correction 5M sl : 

2 

rvTf) 

5MSl = 1R 2 * a (-Ws)ifa bc T^u(p 2 )(f!\ + tfpiMPi)/* (5) 



The s-channel Higgs exchange of Fig. Id gives SM s2,t 

l t I ■ 2 \ ^ob 



am? , . 0x 5ah , ^ 1 



LWim w s w z i=H o h s m i + ""j 1 i (g) 

U(p2)(ffg^u + ffp3^P4u)v(pi) 

The top quark self-energy 5M sel ^ 1 of Fig.le is: 



= Q™t - (-^ s 2 )(T 6 r a ) i ,e(ff 4 )^(p 3 )^(p 2 ) 7M (^2 - h + m t ) 
ISnm^s^ 

[f self,t + _ ^ )]( ^ _ ^ + m<)7 ^ (pi)/( £ _ m? ) 

8M vl,t and SM v2,t of the vertex corrections (Fig.lf-Fig.lg) are: 
+/4 1, *Ta^4 + /s 1 ' t P2 l xh){h + rn t )luv{pi)/{i - m 2 ) 



(7) 



(10) 
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fy 777 

<^ 2 >< = ° « 2 (-^ 2 )(T 6 ni^(P4)"6(p 3 )^(p 2 ) 7 ,(^ - ^ 4 + m t ) 
(A^'Si, + fl U Piu + ff'%7* + ff^faPiMpi)/ (i - mf) 

SM box,t of the box diagram (Fig.lh) is: 

5M b ^ = arr f 2 (-ig s 2 )(^ tt )i^(P4)^(p 3 )^(p 2 )[/ 1 6 ^ + #7,7, 

■ /?6,t -fM x^i* x^M 

+/3 Plulu + Ji VXpTiu + h V2vl» + U V2i.lv + f 7 VlpVlv 

+ft t Pl^P2u + t/V2pV\v + floP2vJ>2u + MflllvlfX- 

i £b,t , rb,t . rb,t rb,t rb,t 

+ Jl2 lulv + JlhPlvln + mVlyTfv + JlkPlvln + JlhP2lJ.lv 

+fviPlnPlv + fliPlnP2v + fl9P2»Plv + fwP2^P2u)]v{pi) 

The following parameters are used in our calculation: 
v/i = UTeV, m t = l76GeV, 

MRS parton distribution set A' [9], Q 2 = s, (11) 
\r]\ < 2.5, p T > 20 GeV. 

The f] and pt cuts are used in order to avoid the singularities in the calculation of 

form factors and increase the relative corrections. 

III. Numerical Results and Discussions 

In our calculation, we checked our program with gauge invariance and it is 
accurate to 10~ 10 . For comparison, we first give the results of the SM in Fig.2. We 
see that the SM results are always negative and do not exceed —1.5%. The shape 
of the curve is like that of Ref. [|J. In Fig. 3 and Fig. 4 we present the results of 
the general 2HDM and the MSSM with a small (3 = 0.25 as in Ref. @. In the 
general 2HDM, the correction can be positive or negative depending on the value 
of triH = TUh- The positive correction can reach about 5% and the negative about 



—7% which may be potentially observable experimentally. The MSSM results are 

all negative and can reach as large as —10% for m H ± < 400 GeV. The results of 

the MSSM with tan (3 = 1 are given in Fig. 5. The corrections never exceed —2%. 

Therefore, if tan f3 cannot be small as may be in the case of the MSSM, the Yukawa 

corrections will be too small to be observable. Although there is a factor — \ tan 2 (5 

mi 

in i] H +, the results vary insensitively with (5 in the large tan/5 region as shown in 
Fig. 6 for the MSSM with tan/? = 70. This conclusion is also valid for a general 
2HDM. 

Therefore, we come to the conclusion that for small /3, in the general 2HDM 
and the MSSM, the Yukawa corrections to gg — > tt at LHC can exceed about —10% 
which may be potentially observable. For larger tan/5 = 1, the Yukawa corrections 
in the general 2HDM and the MSSM will not be observable. 

This work is supported in part by the National Natural Science Foundation of 
China, the Fundamental Research Foundation of Tsinghua University and a grant 
from the State Commission of Science and Technology of China. 
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Appendix 



We give here the form factors for the matrix element appeared in the text. They 
are written in terms of the conventional one-, two-, three- and four-point scalar loop 



integrals defined in Ref . [1(J . We set Mz and Mw to zero in the numerical calculation 
as in Ref. 131. 



A sl = E ^[--2C 24 -^ + m 4 2 (4C -C 21 ) 

i=H0,h A 

-k 2 (C 22 - c 2 3))(-P2,k,m i ,m t ,m t ) 

- E n[l ~ 2C 2 4 - z% - m 2 C 21 

i=A,Z Z 

-k 2 (C 2 2 - C 23 )](-p 2 , k,mi,m t ,m t ) 

E v4-2C 24 -Zv + m 2 (C -C 21 ) 

i=H+,W+ A 

-k 2 (C 22 - C 23 )](-p 2 ,k,mi,m b ,m b ) 
/I 1 = - E Vi[^ m t(2C 11 + C 2 i)] (-p 2 ,k,rrii, mt,m t ) 

i=H°,h 

- E Vi[ 2r ntC 21 ](-p 2 ,k,mi,m t ,m t ) 

i=A,Z 

- E Vi[ 2m t(Cn + C 21 )](-p 2 , k,mi,m b ,m b ) 

i=H+,W+ 

fi 2,t = 4:m t l^ + m 2 Co-p 3 -p4(2C 22 -2C 23 + Co)}(p4;,-k,m t ,m t ,m t ) 
ff' f = 4m t [C + 4(C 22 - C 23 )](p4, ~k, m t , m t , m t ) 
jseif,t _ I ^ r]i[m t (-B + Z- + 5mi)](i,m t ,mi) 

[i=H°,h 

E r]i[m t {B + Z? + 6mi)](i, m t , rrn) 

i=A,Z 

E Vi[ m t(Zi + 8mi)](i,m b ,mi) > /(i-m 2 ) 

i=H+,W+ J 
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ft U = E V l [Bi-Zn(i,m u m, i ) 

[i=H°,h 

E Vi[ B i ~ Z-](i,m t ,mi) 

i=A,Z 

E Vi[Bi - Z?](i,m b ,mi) \ /(t-m 2 t ) 

i=H+,W+ J 

tt U = - E r ll [\-2C 2A -Z v l -m 2 t {2C 11 + C 21 ) 

i=H°,h 

+2p 2 -Pi{Ci2 + C 2 3)](-P2,P4,m i ,m t ,m t ) 
~ E m\\-2C 2A -ZV l -m 2 t {2C ll + C 2l ) 

i=A,Z 1 

+2p 2 -P4(C 12 + C 23 )](-P2,P4, m;,m t ,m t ) 

E ^[^-2C 2 4-^-^t(2Cii + C 21 + C ) 

i=H+,W+ 1 

+2p 2 -Pa{C 12 + C 23 )](-p 2 ,p4,mi,m b ,m b ) 
f 2 l,t = ~ E ?7i[2m t (3Cii + C 2 i + 2C )](-p 2 ,P4,^i,^t,^t) 

- E Vi[2m(C n + C 21 )}(-p 2 ,p A ,mi,m t ,m t ) 

i=A,Z 

- E Vi[2m t (2C u + C 2 i + C )](-p 2 ,p A ,mi,m b ,m b ) 

i=H+,W+ 

jvi,t = ^ 77j[m t (2Co + Cii)](-p 2 ,p 4 ,"7, i ,TO t ,mt) 
E ^[ m <C , ii](-P2,P4,m i ,mi,m t ) 
E Vi[ m t(Cn + C )](-p 2 ,p 4 ,mi,m b ,m b ) 

i=H+,W+ 

fit* = E ^[2(^12 + C23)](-P2,P4,Wj,m t ,m t ) 
E Vi[2(C 12 + C 23 )](-p 2 ,p4,m;,m t ,m t ) 
E Vi[2(C 12 + C 23 )](-P2,P4,mi,m 6 ,m 6 ) 

/r 2 '* = E ^[i-2C 24 -^-m, 2 (2Cn + C 21 ) 

+2pi • J5 3 (C 12 + C 23 )](pi, -p 3 ,m h m t ,m t ) 
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- ]T r ]i [\-2C 2i -Z v i -m 2 {2C 11 + C 2l ) 

i=A,Z 1 

+2pi -p 3 (Ci 2 + C 23 )](pi, -P3,mi,m t ,m t ) 

- r?4-2C 24 -^-™'(2Cii + C 21 + C ) 

i=H+,W+ 

+2pi -p 3 (Ci 2 + C 23 )](pi, -p3,mi,m b ,m b ) 
/2 2 '* = Vil 2m t(3C u + C 2 i + 2C )](pi, -p 3 , rrii, m t , m t ) 

i=H°,h 

Vi[ 2m t(C n + C 2 i)](pi, -p 3: m i: m t: m t ) 

i=A,Z 

Y r)i[2m t (2C u + C 21 + Co)](pi, -p3,rrii,m b ,m b ) 

i=H+,W+ 

jv2,t = _ j2 rji[m t (2C + Cn)](pi, -p 3 , m,, m t , m t ) 

i=H°,h 

~ Yl Vi[ m tC n ](Pi,-P3,mi,m t ,m t ) 

i=A,Z 

- Y Vi[ m t(C 11 + C )](pi, -p3,mi,m b ,m b ) 

i=H+,W+ 

fl 2,t = H ^i[2(Ci2 + C 23 )](pi, -p 3 ,mi,m t ,mt) 

i=H°,h 

r]i[2(C 12 + C 23 )](pi, -p3,m h m t) m t ) 

i=A,Z 

Y Vil 2 (C 12 + C 23 )](pi, -p3,mi,m b ,m b ) 

i=H+,W+ 



f\ l = ~ S ^[/i 2?3 + 4miZ) 27 ](-p 2 ,p4,p 3 ,m i ,mi,m t ,mi) 
- S Vi[ff q \-P2,P±,P3,m,rn u m u m t ) 

i=A,Z 

~ Vi[ff 93 + 2m t D 27 ](-p 2 ,p 4 ,p3,mi,m b ,m b ,m b ) 

i=H+,W+ 

it* = - £ Vilf^ 3 + Zm t {-m 2 D 21 + 2p 2 • p 4 £>24 
+2p 2 • p 3 -D 25 - 2p 3 • p 4 Ase) - 8m t D 27 

+4m t (-m 2 D u + P2 • Pa{Du + D 12 - £>i 3 ))](-p 2 ,P4,p 3 , m*, m t , m t , m t ) 
~ Vilff 93 + m t (-m 2 D 21 + 2p 2 -p 4 D 2i 

i=A,Z 

+2p 2 -p 3 D 25 -2p3 -PaD 2& ) - 4m t D 27 ](-p 2 ,p 4 ,p3,mi,m t ,m t ,m t ) 
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- E Vilff 9 ' + 2m t (-m 2 t D 21 + 2p 2 ■ p 4 D 24 

i=H+,W+ 

+2p 2 ■ p 3 D 25 - 2P3 ■ PiD 26 ) - 6m t D 27 + m t (-3m 2 D u 
+2p 2 -p 4 (D u + Dm - D 13 )) -m 3 t D ](-p 2 ,p 4 ,p 3 ,mi,m b ,m b ,m b ) 
fV = - J2 Vi[ff q3 ](-P2,P4,P3,rni,m t ,m t ,m t ) 

i=H°,h 

~ S Vi[ff 93 }(-P2,P4,P3,mi,m t ,mt,m t ) 

i=A,Z 

- J2 Vilff 93 + 2m 2 t (D 12 - D 13 )](-p 2 ,p 4 ,p 3 ,mi,m b ,m b ,m b ) 

i=H+,W+ 

ff = ~ J2 Vilff 93 -8m 2 D 13 ](-p 2 ,p 4 ,p 3 ,mi,mt,m t ,m t ) 

i=HO,h 

~ J2 Vilff 93 }(-P2,P4,P3,mi,m t ,m t ,m t ) 

i=A,Z 

~ S Vilff 93 - 2m 2 D 13 ](-p 2 ,p 4 ,p 3 ,mi,m b ,m b ,m b ) 

i=H+,W+ 

ff = - J2 Vilff 93 ]{-P2,PA,P^m h m t ,m u m t ) 

i=H°,h 

~ Vilff 93 }(-P2,P4,P3,mi,m t ,m t ,m t ) 

i=A,Z 

~ E Vilff 93 + 2rn 2 (D l2 - D n )](-p 2 ,p 4 ,p 3 ,mi,m b ,m b ,m b ) 

i=H+,W+ 

ff = ~ J2 Vilff 93 + 8m 2 (D + D n - D 13 )](-p 2 ,p 4 ,p 3 ,m i ,m t ,m t ,m t ) 

i=H°,h 

~ E Vilff 93 ](-P2,P4,P3,mi,m t ,m t ,m t ) 

i=A,Z 

~ E Vilff 93 + 2mt(A) + 2^11 - 2D 13 )](-p 2 ,p 4 ,p 3 ,m h m b ,m b ,m b ) 

i=H+,W+ 

fV = ~ E Vilff 93 + 8mtD 2 Q](-p 2 ,p 4 ,p 3 ,m i ,mt,m t ,m t ) 

i=H°,h 

~ Vilff 93 }(-P2,P4,P3,mi,m t ,m t ,m t ) 

i=A,Z 

~ J2 Vilff 293 + 4m t D 26 ](-p 2 ,p 4 ,p 3 ,m i ,m b ,m b ,m b ) 

i=H+,W+ 

ff = ~ Yl Vilff 93 + 8m t (D 26 - D 25 )](-p 2 ,p 4 ,p 3 ,mi,m t ,m t ,m t ) 

i=H°,h 

~ E Vilff 93 ](-P2,P4,P3,mi,m t ,m t ,m t ) 

i=A,Z 

~ Vilff 93 + ^m t (D 26 - D 25 )](-p 2 ,p 4 ,p 3 , m u m b ,m b ,m b ) 

i=H+,W+ 
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/£'* = _ E Vilfg 293 + %rn t {D 2 & - D 24 - D 12 )](-p 2 ,p 4 ,p 3 ,mi,m t ,m t ,m t ) 

i=H°,h 

~ E Vi[ff q3 ](-P2,P4,P3,mi,m t ,m t ,m t ) 

i=A,Z 

~ E Vi[ff q3 + ^ m t(D 26 - D 24 - D 12 )](-p 2 ,p 4 ,p 3 ,mi,m b ,m b ,m b ) 

i=H+,W+ 

fib' = - E Viif!? 3 + 8m t (D 21 + D 26 

i=H°,h 

-D 2A - D 25 + D u - D 12 )] (-p2, Pa, P3, m, m t , m t , m t ) 
~ E Vi[fio q3 K-P2,P4,P3,m,rn t ,rn t ,rn t ) 

i=A,Z 

~ E ^[fw 3 + ^rn t (D 21 + D 26 

i=H+,W+ 

-D 24 - D 25 + D u - D 12 )\ (-p 2 , p 4 , p 3 , rrii, m b , m b , m b ) 
fit = - E Vi[fii q3 ](-P2,P4,P3,rn i ,rn t ,m t ,m t ) 

i=H°,h 

~ E Vi[fii q3 }(-P2,P4,P3,mi,m t ,m t ,m t ) 

i=A,Z 

~ E Vi[fn q3 }(-P2,P4,P3,mi,m b ,m b ,m b ) 

i=H+,W+ 

fi2 = - E ViifiZ 93 + ^m 2 t D ](-p 2 ,p 4 ,p 3 ,mi,m t ,m t ,m t ) 

i=H°,h 

~ E Vi[fi2 q3 ](-P2,P4,P3,m i ,m t ,m t ,m t ) 

i=A,Z 

~ E Vilfu 93 + m 2 t (D + D u - D 13 )] (-p 2 ,P4,P3,mi,m b ,m b ,m b ) 

i=H+,W+ 

fit = - E Vi[fi3 q3 ~ 4m t D 12 ] (-p 2 ,p 4 , Ps, m, m t , m t , m t ) 

i=H°,h 

~ E Vi[fi3 q3 }(-P2,P4,P3,mi,m t ,m t ,m t ) 

i=A,Z 

~ E vAfi? 3 - 2m t D 12 ](-p 2 ,p 4 ,p 3 ,mi,m b ,m b ,m b ) 

i=H+,W+ 

fit = - E Vi[fi4 q3 + ^mD 13 ] {-p 2 ,P4, P3, m h m t , m u m t ) 

i=H°,h 

~ E Vi[f q 4 q3 }(-P2,P4,P3,mi,m t ,m t ,m t ) 

i=A,Z 

~ E Vi[f i4 93 + 2m t D 13 ](-p 2 ,p 4 ,p 3 ,mi,m b ,m b ,m b ) 

i=H+,W+ 

fit = - E Vi[f i5 q3 + ^ m t(D n ~ D 12 )](-p 2 ,p 4 ,p 3 ,mi,m t ,m t ,m t ) 

i=H°,h 
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- Z! Vi[fi5 q3 }(-P2,P4,P3,mi,m t ,m t ,m t ) 

i=A,Z 

~ J2 Vi[f?5 q3 + 2m t (D n - D 12 )](-p 2 ,P4,P3,m,m b ,m b ,m b ) 

i=H+,W+ 

fw = - J2 Vi[fi6 qS + 4m t (Di 3 - D 11 )\(-p 2 ,p i ,p 3 ,m i ,rn t ,m u rn t ) 

i=H°,h 

~ Z Vi[fi6 q3 ](-P2,P4,P3,m,rn t ,m t ,m t ) 

i=A,Z 

~ Z Vilfw 93 + 2m t(D 13 - Du)](-p2,P4,P3,rni,m b ,m b ,m b ) 

i=H+,W+ 

fi7 = ~ J2 Vi[f q 7 q3 ](-P2,P4,P3,mi,m t ,m t ,mt) 

i=H°,h 

~ J2 Vi[fi7 q3 ](-P2,P4,P3,mi,m t ,m t ,m t ) 

i=A,Z 

~ Vi[fi7 93 ](-P2,P4,P3,mi,m b ,m b ,m b ) 

i=H+,W+ 

fit = ~ 2 Vi[f q 8 q3 K-P2,P4,P3,m,rn t ,rn t ,rn t ) 

i=H°,h 

~ 2 Vi[fi8 q3 K-P2,P4,P3,m,rn t ,m t ,m t ) 

i=A,Z 

~ Vi[fi8 q3 }(-P2,P4,P3,mi,m b ,m b ,m b ) 

i=H+,W+ 

fig = ~ J2 Vi[f q 9 q3 ](-P2,P4,P3,rn i ,rn t ,m t ,m t ) 

i=H°,h 

~ Vi[fi9 q3 K-P2,P4,P3,rn i ,rn t ,m t ,m t ) 

i=A,Z 

~ Z Vi[fi9 93 }(-P2,P4,P3,m i ,m b ,m b ,m b ) 

i=H+,W+ 

fw = ~ Z Vi[f20 93 ](-P2,P4,P3,m,rn t ,rn t ,m t ) 

i=H°,h 

- Z Vi[f q o q3 ](-P2,P4,P3,m,rn t ,m t ,m t ) 

i=A,Z 

- Yl Vi[f2o q3 ](-P2,P4,P3,mi,m b ,m b ,m b ) 

i=H+,W+ 

3 = 2m t D 311 

jq2q3 = _ m 3 + 2m t (p 2 • P4-D34 + £>2 ' ^3-^35 ~ P3 " P4D3IO ~ 2-D 3 ll) 

-m 3 t D 21 + 2m t [p 2 ■ p 4 (D 21 - D 25 ) - p 2 ■ p 3 D 25 + p 3 ■ p 4 D 26 + D 27 ] 
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ff g3 = 2[m 2 (D 34 -D 35 ) + 2p 2 -p 3 {D 37 -D 3W ) + 2p 3 -p 4 (D 38 -D 39 ) 

+2p 2 ■ p 4 (D 310 -D 36 ) + 4D 312 - 6D 313 ] 

+4m 2 (D 24 - D 25 ) + Ap 2 ■ p 4 (D 25 + D 26 - D 22 - D 23 ) + 4D 27 
ff q3 = 2[m 2 D 35 -2p 2 -p 3 D 37 + 2p 3 -p 4 D 39 -2p 2 -p 4 D 310 + 4D 313 \ 

+%> 2 -Pa{D 23 - D 2b ) 
ff q3 = 2[m 2 (D 34 -D 31 ) + 2p 2 -p 4 (D 34 -D 36 )+2p 2 -p 3 (D 35 -D 310 ) 

+2p 3 ■ p 4 (D 38 - D 310 ) - 4D 3n + 4D 312 ] 

+4m 2 (D 24 -D 21 )+ 4p 2 • p 4 (D 24 - ^22) 
ff q3 = 4(^311 - £313) + 2m 2 (D 21 - 2D 25 ) + 4p 2 ■ p 3 D 25 + 4p 2 ■ p 4 D 26 - Ap 3 ■ p 4 D 26 
ff q3 = 4m t D 3W 
ft 3 = 4m t (D 310 -D 35 ) 
ff q3 = 4m t (D 310 -D 34 )-4m t D 24 
f q2q3 = 4m t (D 31 -D 34 -D 35 + D 310 )+4m t (D 21 -D 24 ) 
f q2q3 = 2(D 313 -D 312 )-2D 27 

fu 3 = rn 2 (D 34 - D 35 ) + 2[p 2 ■ p 4 (D 310 - D 36 ) + p 2 ■ p 3 {D 37 - D 310 ) 
+Ps ■ Pa(D 38 ~D 39 ) + 2D 312 - 2D 313 ] 
+m 2 (2D 24 - 2D 25 - D 21 ) 

+2[p 2 ■ p 4 (D 25 + 2D 26 - D 22 - D 23 ) + p 2 ■ p 3 D 25 - p 3 ■ p 4 D 26 ] 
f q l q3 = -2m t D 24 
ff 4 q3 = 2m t D 25 
f q2q3 = 2m t (D 21 -D 24 ) 
f q2q3 = 2m t (D 25 -D 21 ) 
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= 4(^39 -^38)+ 4(^23-^26) 

fis q3 = 4(D 39 + D 3W -D 37 -D 38 )+A(D 25 -D 26 ) 

/g ?3 = 4(^36 + ^39 " ^38 - D 3W ) + 4(^22 + D 23 - D 25 - D 26 ) 

/f 2<z3 = A(D 35 + D 36 -D 34 -D 37 -D 38 + D 39 )+A(D 22 -D 24 + D 25 -D 2e ). 



In the above formulae 

k = pi +P2 = P3 +P4, s = k 2 , i=(p 2 -p A ) 2 , u = (p 2 -p 3 ) 2 , 

T» = (- P4+P3 Ye(p 3 ) • e(p 4 ) + (2p 4 +Ps) ■ e{p 3 )e{p A Y ~ (2p 3 + Pa) • e(p 4 )e(p 3 ) M , 
sin 2 ct cos 2 a 

sin p sin 5 

2 

?7 A = cot 2 /3, 77^+ = (cot 2 j3 H |tan 2 /3), r?z = W+ = 1- 

The renormalization constants are: 



Z\ = [B Y + 2m 2 (G + GOKm^m^mO /or i = H°, h 
Zl = [B, + 2m 2 (G 1 - G )](m 2 ,m t ,m t ) for 1 = A, Z 
Z\ = [B, + 2m 2 G 1 ](m 2 ,m b ,m l ) for % = H+, W + 

Srrii = [Bq — Bi](m 2 ,mt,mi) for i = H°, h 
Srrii = [-B Q - B^{m 2 t ,m u m^) for i = A, Z 
5rrii = [— i?i](m 2 , rrib, m^) for i = H + , W + 

, n dB (p 2 ,m 1 ,m 2 ) _ dB 1 (p 2 ,m 1 ,m 2 ) 
wnere Lr — ^ | p 2 =m 2, — — | p a =m 2 
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Figure Captions 

Fig.l Feynaman diagrams of tree level and 0{am 2 t jm^ v ) Yukawa corrections of gg — > 
tt. 

Fig. 2 Results of Yukawa corrections in the Standard Model. 

Fig. 3 Results of Yukawa corrections in the general Two-Higgs-Doublet Model with 
P = 0.25. 

Fig. 4 Results of Yukawa corrections in the Minimal Supersymmetric Standard Model 
with = 0.25. 

Fig. 5 Results of Yukawa corrections in the Minimal Supersymmetric Standard Model 
with p — 7r/4. 

Fig. 6 Results of Yukawa corrections in the Minimal Supersymmetric Standard Model 
with tan/3 = 70. 
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